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Abstract. We report ISO-LWS far infrared observations of CO, water and oxygen lines towards the protobinary 
system IRAS4 in the NGC1333 cloud. We detected several water, OH, CO rotational lines, and two [OI] and 
[CII] fine structure lines. Given the relatively poor spectral and spatial resolution of these observations, assessing 
the origin of the observed emission is not straightforward. In this paper, we focus on the water line emission 
and explore the hypothesis that it originates in the envelopes that surround the two protostars, IRAS4 A and 
B, thanks to an accurate model. The model reproduces quite well the observed water line fluxes, predicting a 
density profile, mass accretion rate, central mass, and water abundance profile in agreement with previous works. 
We hence conclude that the emission from the envelopes is a viable explanation for the observed water emission, 
although we cannot totally rule out the alternative that the observed water emission originates in the outfiow. 
The envelopes are formed by a static envelope where the density follows the law, at r > 1500 AU, and a 
collapsing envelope where the density follows the r~^^'^ law. The density of the envelopes at 1500 AU from the 
center is ~ 4 x 10^ cm~^ and the dust temperature is ~ 30 K, i.e. about the evaporation temperature of CO-rich 
ices. This may explain previous observations that claimed a factor of 10 depletion of CO in IRAS4, as those 
observations probe the outer < 30 K region of the envelope. The water is ~ 5 x 10~^ less abundant than H2 in the 
outer and cold envelope, whereas its abundance jumps to ~ 5 x 10~® in the innermost warm region, at r < 80 AU 
where the dust temperature exceeds 100 K, the evaporation temperature of H20-rich ices. We derive a mass of 0.5 
Mq for each protostar, and an accretion rate of 5 x 10~^ Mq yr~^ , implying an age of about 10000 years, if the 
accretion rate remains constant. We finally discuss the difference between IRAS4 and IRAS16293-2422, where a 
similar analysis has been carried out. We found that IRAS4 is probably a younger system than IRAS16293-2422. 
This fact, coupled with the larger distance of IRAS4 from the Sun, fully explains the apparent difference in the 
molecular emission of these two sources, which is much richer in IRAS16293-2422. 
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1. Introduction 

The south part of the NGC1333 reflection nebulae, in the 
Perseus cloud, is an active star forming region, contain- 
ing many infrared sources associated with molecular flows 
and numerous Herbig-Haro objects. IRAS4 was flrst iden- 
tified by Jennings et alj ( |1987 ), and further observations 
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(Sandell et al. 1991) revealed IRAS4 it was a binary sys- 
tem resolved into two components, named IRAS4A and 
IRAS4B, and separated by 31". Interferometric observa- 
tions (Lay et al. 1995; Looney et al. 2000) have shown 



further multiplicity of the two sources. IRAS4A is itself a 
binary system with a separation of 10", and there is some 
evidences that IRAS4B could also be a multiple system, 
with a separation of 0.5". 

The distance of the NGC1333 cloud is much debated. 



Herbig & Jones (1983) found a distance of 350 pc for the 



Perseus 0B2 association (a more recent estimate based 
on the Hipparcos data gives 318±27; ie Zeeuw et al 



1999|) , but ext inction observations towards NGC1333 itself 



( |Cernis||l990[ ) suggest that it may be as clos e as 2 20 pc. 
Assuming a distance of 350 pc, ^andell et al. (1991) mea- 
sured a system total luminosity of 28 Lq (11 Lq at 220 
pc) equally shared between IRAS4A and B. They derived 
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an envelope mass of 9 and 4 Mq respectively (3.5 and 
1.5 Mq at 220 pc). This relatively large mass, together 
with the low bolometric luminosity suggest that both 
sources are deeply embedded and probably very young. 



They have been classified as Class sources (Andre et al 



1993). IRAS4A and B are both associated with molecu- 



lar outflows, detected in CO, CS (Blake et al. 1995) and 
SiO (Lefloch et al, 1998) millimeter transitions. The out- 
flow originating from IRAS4A is very highly coUimated, 
whereas that originating from IRAS4B is rather compact 



and u nresolved in single dish observations ( Knee & Sandell 



2000| ). The dynamical ages of both outflows are a few thou- 
sands years. 

In the past years, many observational studies have 
been focused on the continuum emission of IRAS4. Recent 



et al. 1996, hereafter CHT96; Doty & Neufeld 1997). Given 
the relatively large range of level energies (from ^ 100 
to ^ 500 K) and spontaneous emission coefficients (from 
10^^ to ~ 1 s~^) of the water transitions observed by 
ISO-LWS, the observed lines can in particular probe the 
innermost regions of the envelope. This makes the analysis 
of the ISO-LWS water lines a precious and almost unique 
tool (when considering the water abundance across the en- 
velope). The reverse of the coin is that assessing the actual 
origin of the water emission is somewhat difficult and still 
debated, as the spectral and spatial resolutions of ISO- 
LWS are relatively poor to disentangle the various com- 
ponents falling into the beam. For example, strong molec- 
ular line emission is often associated with the outflows 



emanating from young protostars (e.g. Bachillcr fc Perez 



works nclude maps of the region obtained with IRAM at Gutierrez 1997 ). As already mentioned, the line emission 



1.3 mm ( Lefloch et al.||l998|) an d with SCUBA at 450 and 



850 fiia (Sandell fc Knee 2001). An accurate modeling of 



the continuum emission has been very recently carried out 
by l0rgensen et al. ( [2002 , hereafter JSD02), who recon- 
structed the dust temperature and density profiles across 



the two envelopes. 



The molecular line emission is probably a better and 
certainly a complementary tool to probe the dynami- 
cal, chemical and physical structure of the envelopes of 
IRAS4. The last decade has seen flourishing several stud- 
ies of molecular line profiles (e.g. Gregersen et al. 1997; 
Evans 1999) and line spectra (Blake et al, 1995D , all hav- 
ing in common the goal of reconstructing the physical 



from CO, CS and other molecules are certainly contam- 
inated by the outflowing gas in IRAS4. Nonetheless, low 
lying lines seem to be more affected than high lying lines 
in first instance, and different molecules suffer differently 
from this "contamination" , as proved by the Di Francesco 
et al.| ( |200l| ) observations. Although water has been pre- 



dicted to be very abundant in shocked gas, the published 
ISO observations show that the water emission is usu- 
ally stronger towards the central sources and weaker, if 
detected at all, in the direction of the peaks of the out- 



flows powered by low mass protostars (see Ceccarelli et al 
pOOOa , for a review 



structure of the protostellar envelopes. Specifically, Blake 



When water lines are detected in 
clear-cut shocked regions, the water abundance seems to 
be lower than that predicted by the models, like in the 



et al. ( L995) carried out a multifrcquency study of several 



moleci les in IRAS4, including H2CO and CH3OH. Their 



two major results regarding the structure of the IRAS4 en- 
velopes are: 1) a large depletion, around a factor 10-20, of 
CO and all molecules in the envelope, and 2) the presence 
of a region with an increased abundance of CS, SiO and 
CH3OH, that the authors attribute to mantles desorption 
caused by grain-grain collisions induced by the outfiows 
originating from the two p rotostars. More rec e ntly i nter- 
fero metric observatio ns by Di Francesco et al. (2001) (see 
also Choi et al. 1999) detected an inverse P-cygni profile 
of the H2CO 32.1 — 2i_i line on a 2" scale towards both 
IRAS4A and B, providing the least ambiguous evidence of 
infall motion towards a protostar ever. From a simple two- 
layer modeling, they derived an accretion rate of 1.1 x 10^"* 
and 3.7 x 10~^ Mq yr~^ , an inner mass of 0.7 and 0.2 
Mq , and an age of 6500 and 6200 yr (assuming constant 
accretion rate) for IRAS4A and IRAS4B respectively. 

In this paper we concentrate on the far infrared (FIR) 
line spectrum, and in particular the water line spectrum 
observed with the Long Wavelength Spectro meter ([Clegg 



case of HH54 ([Liseau et al\ |1996| ) or HH7-11 ( [Vlolinar: 
et al.||199g| , |2000D , or in the outfiows of IRAS4 (see next 



et al. 1 396, herein after LWS) on board ISO (Kessler et al 



1996) in the direction of IRAS4. The goal of this study is 



to check whether the observed water line emission can be 
attributed to the thermal emission of the envelopes sur- 
rounding the IRAS4 protobinary system. Water lines have 
in fact been predicted to be a major coolant of the gas in 
the collapsing envelopes of low-mass protostars (Ceccarelli 



section). Finally, SWAS observations seem to support the 
evidence that the w ater abundance in th e shocked regions 
is a few times 10^^ ( Neufeld et al. 200C). These facts, to- 
gether with the apparent correlation between the observed 
water emission and the 1.3 mm continuum, and the lack 
of correlation with SiO emission]^ in low mass protostars 
( Ceccarelli et al. 1999| ) play in favor of a relatively low 
contamination of the ISO-LWS observed water emission 
by the outflow and encourage us to explore in detail this 
hypothesis for the IRAS4. 

In the specific case of IRAS4, the Suhmillimiter 
Wavelengths Astronomical Satellite (SWAS; Melnick et al, 
^000| ) observed the ground 0-H2O line at 557 GHz 
(Neufeld et al. 2000, Bergin et al. 2002). Given its rel- 
atively large linewidth ('^18 km s^^) the 557 GHz line is 
certainly dominated by the outfiow emission. Nonetheless, 
this does not imply that the ISO FIR water lines also orig- 
inate in the outflow, and this for two reasons. First, the 
beamwidth 4') of the SWAS observations, being about 
3 times that of ISO-LWS, encompasses the entire outflow, 
whereas the ISO observations do not encompass the two 



SiO emis sion is usually associated with the outflow strong 



shocks (e.g. BachiUer fc Perez Gutierrez 1997), and it is be- 



lieved to be a product of grain mantle desorption ( Caselli et al 
|l997| ; [Schilke et aL||l997| ) 
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emission peaks of the outflow (see also 3.1), but only the 
envelope. Second, the 557 GHz transition, being the wa- 
ter ground transition, is more easily excited than the FIR 
water lines, and therefore the latter probably probe dif- 
ferent regions. In fact, Bergin et al. (2002) find that most 
of the 557 GHz line must originate in a component colder, 
hence different, by that probed by the FIR water lines, 
even under the assumption that they probe the outflow. 
To summarize, decicing whether the observed FIR water 
emission in IRAS4 originates in the outflow or in the en- 
velope remains an open question, based on the available 
present observations. In this article we explore in detail 
the latter hypothesis and submit it to the scrutiny of an 
accurate modeling, trying hence to answer to the ques- 
tion on a theoretical basis. At this scope we used the 
CHT96 model, already successfully applied to the solar 
type protostar IRAS16293-2422 which allowed to explain 
more than two dozen observed ISO-LWS water lines and 



grou nd-base d millimeter SiO and H2CO lines (Ceccarelli 
et al. | B0004 |b|)- One of the major results of that work is 



the prediction of the existence of a hot core like region in 
the innermost part of the envelope of IRAS16293-2422, in 
which the dust temperature exceeds the evaporation tem- 
perature of interstellar ice (~ 100 K). These studies have 



been confirmed by the recent analysis by Schoier et al 
( 2002| ) of several other molecular transitions. Such hot 



cores are well studied around massive protostars where 
- driven by reactions among the evaporated ice molecules 
in the warm gas - their chemical composition differs sub- 



stantially from that of quiescent clouds (Walmsley 1989; 



Charnley et al. 1992). Hot cores around low mass proto- 



s tars may actual l y have a different chemical composition 
( [Ceccarelli et al. 2000b). This molecular complexity may 
be of prime interest on account of a possible link to the 
chemical history of the solar nebula and hence the molec- 
ular inventory available to the forming Earth and other 
solar system planets and satellites. 

In order to understand the physical and chemical pro- 
cesses that take place during the first stages of star for- 
mation, it would be necessary to undertake a work similar 
to that the one done on IRAS16293-2422 on a larger sam- 
ple of protostars. In this paper we present a study of the 
structure of the envelope of NGC1333-IRAS4, obtained 
using ISO-LWS observations of the H2O far-infrared lines. 
A preliminary analysis of the same set of data has al- 



ready been presented in Ceccarelli et al. (1999) and Caux 
et al. | (|l999aD . Here we revisit the data using a new cal- 



ibration and compare the observations with the CHT96 
model predictions, testing a large range of model parame- 
ters. This study is part of a large project aimed to model 
the water emission in several low mass protostars. The wa- 
ter observations are complemented with formaldehyde and 
methanol ground based observations, to have a complete 
budget of the most abundant molecules in the innermost 
regions of the protostellar envelopes (Maret et al. in prepa- 
ration). Finally, the structure obtained by the analysis of 
these observations will be compared with that indepen- 
dently obtained by continuum observations by JSD02. 



The outline of the article is the following. In ^ we 
present the data, in 33 we describe the modeling of the ob- 
served lines and in §Erwe discuss the physical and chemical 
structure of the envelope, namely the density and tem- 
perature profiles, as well as the abundances of the major 
species across the envelope. Besides, the central mass of 
the protostar and its accretion rate can also be constrained 
by these observations and modeling, yielding an alterna- 
tive method to measure these two key parameters. In 
we compare the results of the present study with previ- 
ous studies of IRAS4. Finally, we discuss the similarities 
and differences between IRAS4 and IRAS16293-2422, and 
highlight the importance of complementary ground-based, 
higher spatial and spectral resolution observations to un- 
derstand the physical and chemical processes taking place 
in the innermost regions of low-mass envelopes. 

2. Observations and results 

A full-range spectral survey (43-196 /im) of the IRAS4 
region was performed using LWS. The observations were 
obtained on three positions in LWS grating mode, with 
a spectral resolution of about 200 (AOT LOl). The first 
position was centered in between IRAS4A and IRAS4B 
(02000 = 03^29'"11.9^ (52000 = 31°13'20.3"), and hence 
LWS 80" beam includes both sources. The two other 
positions aimed to the lobe peaks of the outflow pow- 
ered by IRAS4: NE-red, aaooo = 03'^29'"15.6^ (52ooo = 
31°14'40.1"), and SW-blue, azooo = 03'^29"06.6", 
(52000 = 31°12'08. 7"). These observations, performed dur- 
ing revolution 847, are made of 30 scans on the central po- 
sition and 10 scans on the two other positions. The sam- 
pling rate was 1/4 of the grating resolution element (0.29 
/im in the 43-92 fim range, and 0.15 /im in the 84-196 /im 
range). The integration time for each sampled point was 
12 sec on the central position, and 4 sec on the outflow 
positions. 

The same data set has been previously analyzed by 
[Ceccarelli et al. (1999) and Giannini et al. ( 2001 , here- 
after GNLOl). We here report again the analysis of this 
dataset as we used an improved data processing module, 
allowing the correct ion of the t ransient effects that affect 
the LWS detectors ( [Caux[pb0l[ ). As for standard pipeline 
products, the data are calibrated against Uranus, and the 
calibration uncertainty is estimated to be better than 30% 
(Swinyard et al. 1998). On the three observed positions, 
the spectra were then defringed for all ten detectors, and 
the continuum was removed fitting a polynomial baseline 
outside the lines. The data were then averaged over the 
ten detectors and binned at a single resolution to pro- 
duce a single spectrum for each observed position. The 
line fiux measurements were finally performed with the 
ISAP package, using gaussian fits. A particular attention 
was given to the determination of uncertainties associated 
with the line fluxes measurement. These are due to the 
statistical uncertainties, the absolute calibration and the 
baseline determination uncertainties that affect low reso- 
lution LWS observations. In grating mode, on a spectrum 
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rich in lines as it is the case here, the hne confusion may 
be important, making the basehne determination difficult. 
This uncertainty, often neglected in the literature, should 
be taken into account, as it can lead to important errors, 
especially for faint lines. 

Figure |^ shows the observed 60-200 ^m spectra in the 
three observed positions, and Table 1 reports the mea- 
sured line fluxes. The errors quoted in this table include 
statistical errors, errors due to the uncertainty of the base- 
line removal, and an absolute calibration error of 30%. 
The first striking result of these observations is the dra- 
matic difference within the three spectra: while that in- 
cluding IRAS4A and 4B is very rich in CO and water lines, 
molecular emission is barely detected towards the outfiow 
peaks (where the millimeter CO emission is the bright- 
est). On the contrary, the fine structure [01] 63/im and 
[CII] 157/xm lines have comparable fluxes in the three ob- 
served positions. Finally, we wish to comment our results 
with respect to previous published data reductions. The 
present line flux determination agrees with that quoted 
by Ceccarelli et al. (199!;) when considering the uncer- 
tainties. However, we note several differences with respect 
to the values quoted in GNLGl. While there is a relatively 
good agreement between their fluxes of the strongest lines 
and ours (see Table ^ , there is a noticeable discrepancy 
between our respective reductions regarding the weakest 
lines. We think that this is probably due to a too opti- 
mistic evaluation of the noise in GNLGl . For example, in 
the NE-red position we find a statistical error around 175 
fim of 1.5x 10^^'^ erg s^^ cm^^ /im^^, while Giannini et 
al. quote 2 x 10"^** erg cm~^ ij,m~^ in their Table 3. 
We do not confirm neither the detection of CO lines with 
J«p > 2l|, nor the 125.4, 83.3, 66.4 and 58.7 ^m water 
lines on-sourcc. As shown in our Fig. |^, we only detected 
the 179.5 /xm 174.6 /im and 108.0 /im lines in the outflow 
peak position NE-red, and the CO lines between Jup — 14 
and 17. We also do not conflrm their detections of CO J„p 
> 18 lines in the NE-red outflow peak position. Finally, in 
the SW-blue position we only detected C+ and 01 63 /im 
emission and very marginally the H2O 179 /im lines. 



3. Line modeling 

3.1. Origin of the FIR line emission 

The molecular emission (H2O, CO and OH) observed to- 
ward IRAS4 can have at least three different origins: the 
two outflows powered by IRAS4A and IRAS4B, the PDR 
at the surface of the cloud, and the collapsing envelopes 
around the two protostars. The origin of the molecular line 



^ CO 21-20 is definitively undetected, CO 22-21 is blended 
with OH at 119 /xm, while CO 23-22 and CO 24-23 are blended 
with water lines at 113.5 and 108.8 /xm line respectively. One 
can note that the 108 /xm is wider than expected, which could 
be due to the presence of the CO 24-23 line. Nonetheless, this 
detection would be inconsistent with the non detection of lower 
and higher J CO lines. Higher spectral resolution observations 
are needed to settle the question. 



emission can be disentangled when the spatial distribu- 
tion of the line emission and/or the line proflles are avail- 
able. For example, lines arising in the envelope or in the 
molecular cloud have narrow proflles whereas lines aris- 
ing from outflows show broadened proflles with extended 
wings. Unfortunately, in the case of the ISO-LWS obser- 
vations, the relatively low spatial and spectral resolution 
do not allow to observationally disentangle the different 
components. However, the comparison between the central 
and the NE-red and SW-blue positions allows a flrst guess 
of the origin of the observed emission. The [01] and [CII] 
lines have comparable line fluxes in the three observed po- 
sitions. For this reason it is likely that the observed [01] 
and [CII] emission is associated with the ambient diffuse 
gas, either emitted in the PDR or in the molecular cloud 
itself. On the contrary, only the lowest lying {Jup < 17) 
CO lines are detected on the NE-red position, while no 
H2O or CO emission is clearly detected on the SW-blue 
position. On the other hand, the observations of the CO 3- 
2 line show that the high velocity gas (the fastest outflow 
component) peaks at the NE-red and SW-blue positions 
(Blake et al. ^995). The lack of water emission in these 
two outflow peak positions is not in favor of the hypothe- 
sis that the on-source water emission originates in the out- 
flow. Although we cannot exclude a different origin and/or 
contamination for example from the densest parts of the 
outffow located in the ISO beam, in the following we ex- 
plore the hypothesis of the envelope thermal emission and 
interpret the observed water line emission according to the 
CHT96 model. The first goal of our modeling is to verify 
that the thermal emission from the surrounding envelope 
can reproduce the water line observations, a necessary con- 
dition even though not sufficient to test this hypothesis. 
A following section will then address the possible origin of 
the CO and OH observed emission. 



3.2. Model description 

The CHT96 model computes in a self-consistent way the 
radiative transfer, thermal balance, and chemistry of the 
main gas coolants (i.e. O, CO and H2O) across the en- 
velope, in the inside-out framework ( ^hu| 1977). Here we 
give a brief description of the main aspects of the model. 

The initial state of the envelope is assumed to be an 
isothermal sphere in hydrostatic equilibrium, which den- 
sity is given by: 

„2 



"-H2 (r) 



(1) 



27r /imnG 

where a is the sound speed, mn is the hydrogen mass, 
fi is the mean molecular mass, r the distance from the 
center and G the gravitational constant. 

At t = the equilibrium is perturbed and the collapse 
starts from inside out, propagating with the sound speed. 
The density in the inner collapsing region is given by the 
free-fall solution: 



1 



M 



ATTumn (2GM,)i/2 



,,-3/2 



(2) 
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Specie 


Transition 


Wavelength 


E 




Fluxes 








(fim) 


(cm-1) 


On source 


NE-red 


SW-blue 


0-H2O 


221-221 


180.49 


134.9 


1.1 ±0.5 


< 0.5 


< 0.5 




2i2-loi 


179.53 


79.5 


2.7 ± 1.0 


0.5 ±0.4 


0.6 ± 0.4 




3o3-2l2 


174.63 


136.7 


1.6 ± 0.7 


0.5 ± 0.4 


< 0.5 




4l4-3o3 


113.54 


224.5 


3.0 ± 1.0 


< 0.5 


< 0.5 




221-I10 


108.07 


134.9 


2.0 ± 0.8 


0.9 ± 0.5 


< 0.5 




5o5-4l4 


99.48 


325.3 


1.4 ± 0.6 


< 0.5 


< 0.5 




6l6-5o5 


82.03 


447.3 


1.2 ± 0.6 


< 0.5 


< 0.5 




423-3l2 


78.74 


300.5 


1.4 ± 0.6 


< 0.5 


< 0.5 




321-2l2 


75.38 


212.1 


1.7 ± 0.7 


< 0.5 


< 0.5 


P-H2O 


322-3l3 


156.19 


206.3 


0.5 ±0.4 


< 0.5 


< 0.5 




3l3-2o2 


138.53 


142.3 


1.2 ± 0.6 


< 0.5 


< 0.5 




220-111 


100.98 


136.2 


1.0 ± 0.5 


< 0.5 


< 0.5 




322-2ll 


89.99 


206.3 


1.4 ± 0.6 


< 0.5 


< 0.5 




331-220 


67.09 


285.1 


0.6 ± 0.4 


< 0.5 


< 0.5 


CO 


14-13 


186.00 


403.5 


1.3 ± 0.6 


1.2 ± 0.4 


< 0.5 




15-14 


173.63 


461.1 


1.5 ± 0.7 


0.8 ± 0.4 


< 0.5 




16-15 


162.81 


522.5 


1.7 ± 0.7 


0.7 ± 0.4 


< 0.5 




17-16 


153.26 


587.7 


1.7 ± 0.7 


0.6 ± 0.5 


< 0.5 




18-17 


144.78 


656.8 


1.3 ± 0.6 


< 0.5 


< 0.5 




19-18 


137.20 


729.7 


1.3 ± 0.6 


< 0.5 


< 0.5 




20-19 


130.37 


806.4 


1.2 ± 0.6 


< 0.5 


< 0.5 


OH 


^113/2,5/2-^113/2,3/2 


119.33 


83.7 


1.0 ± 0.5 


< 0.5 


< 0.5 




113/2,7/2- 113/2,5/2 


84.51 


201.9 


1.1 ± 0.5 


< 0.5 


< 0.5 




ni/2.1/2- 113/2,3/2 


79.15 


126.3 


1.1 ± 0.5 


< 0.5 


< 0.5 


[OI] 


'Pi-'Po 


145.48 


227.7 


< 0.5 


< 0.5 


< 0.5 




^P2-^P2 


63.17 


158.7 


2.2 ± 0.9 


4.3 ± 1.5 


2.2 ± 0.9 


[CII] 


P3/2- Pl/2 


157.74 


63.7 


2.2 ± 0.9 


1.8 ± 0.7 


1.1 ± 0.5 



Table 1. Measured line fluxes in units of 10 ergs s ^ cm ^. Upper limits are given as 2 a. 



where M, is the star mass, M is the accretion rate, 
related to the sound speed by: 



M = 0.975 



G 



v{r) = 



The free-fall velocity is given by: 



(3) 



(4) 



The gravitational energy is released as material falls at 
the core radius i?* , so that the luminosity of the protostar 
is: 



(5) 



In the following is the bolometric luminosity of 
IRAS4, and we leave M, and M as free parameters. 

The radiative transfer in the envelope is solved in 
the escape probability approximation in presence of warm 
dust, following the Takahashi et al. ( 1983| ) formalism. The 
CHT96 model assumes that the initial chemical composi- 
tion is that of a molecular cloud, and then it solves the 
time dependent equations for the chemical composition of 
44 species, as the collapse evolves. H2O, CO and O are of 
particular importance since they are the main coolants of 



the gas, and hence we study the chemistry of these species 
in detail. The CO molecule is very stable, and its abun- 
dance results constant across the envelope. H2O is mainly 
formed by dissociative recombination of the HaO"*" in the 
cold outer envelope, while, at dust temperature above 100 
K, icy grain mantles evaporate, injecting large amounts 
of water into the gas phase. When the gas temperature 
exceeds ~250 K, the H2O formation is dominated by the 
endothermic reactions O ± H2 ^OH ± H followed by H2 
± OH H2O ± H, which transform all the oxygen not 
locked in CO molecules into H2O. 

From the above equations and comments, the water 
line emission depends on the mass of the central object, 
the accretion rate and the abundance of H2O in the outer 
envelope and in the warm region, where its abundance is 
dominated by the mantle evaporation. All these quanti- 
ties directly enter into the H2O line emission, and specif- 
ically into the determination of the H2O column density. 
In fact, the accretion rate sets the density across the pro- 
tostellar region (Eqs. |^ and The central mass of the 
protostar affects the velocity field, and hence indirectly 
the line opacity (Eq. ^) . This parameter also sets the den- 
sity in the free-fall region (Eq. ^) and therefore the gas 
column density in this region. The water emission also 
depends indirectly on the O and CO abundances, which 
enter in the thermal balance and hence in the gas tem- 
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perature determination. Several recent studies (Baluteau 
et al. 



1997 



Caux et al. 1999b; Vastel et al 



2000 



Lis et al. 



2001) have shown that almost the totality of the oxygen 
in molecular clouds is in atomic form. Accordingly, we 
assume the oxygen abundance to be the standard inter- 
stellar value, i.e. 5 x lO^"' with respect to H2. With regard 
to the CO abundance, following Blake et al, (1995) we 
adopt a CO abundance of 10^'°^ with respect to H2, lower 
than the standard abundance as this molecule is believed 
to be depleted on IRAS4. We will comment later on the 
influence of these parameters on our results. Finally, the 
water abundance in the cold and in the warm parts of the 
envelope are poorly known and are free parameters in our 
study. 

To summarize, we applied the CHT96 model to IRAS4, 
and to reproduce the observations we varied the four fol- 
lowing parameters: the mass of the central object M^, , the 
accretion rate M , the water abundance in the outer cold 
envelope Xgut , and the water abundance in the region 
of mantle evaporation Xin ■ The principal limitation to 
the application of this model to IRAS4 is that the ISO 
beam includes both IRAS 4A and IRAS 4B envelopes. As 
a first approximation, we assumed that the two envelopes 
contribute equally to the molecular emission. Finally we 
assumed that the two envelopes touch each other, namely 
they have a radius of 3000 AU (i.e. 30" in diameter), in 
agreement with millimeter continuum observations ( Mottc 
& Ancr^ 2001). In our computations we adopted a dis- 
tance of 220 pc in agreement with JSD02 and a luminosity 



of 5.5 Lq for each protostar, according to Sandell et al 
( |1991| ) when assuming such a distance. 

3.3. Water line modeling results 

In order to constrain the mass, accretion rate and wa- 
ter abundance across the envelope, we run several models 
varying the central mass from 0.3 to 0.8 Mq , the accretion 
rate from 10~^ to 10"'* Mq yr^^ , the water abundance 
in the outer parts of the envelope Xout between 10~^ and 
10~^, and a water abundance in the inner parts of the 
envelope Xin between 10~^ and 2 x 10"^ respectively. In 
the following we discuss the results of this modeling. 

3.3.1. Accretion rate and water abundance in the 
outer parts of the envelope 

One of the difficulties in constraining the central mass, ac- 
cretion rate and water abundance in the envelope is that 
the water line intensity a priori depends on all the pa- 
rameters. However, choosing appropriate lines can help 
constraining one parameter at once. Low lying H2O lines 
are expected to rapidly become optically thick in the outer 
envelope, where they are easily excited. Hence these lines 
depend weakly on Xin and M, (which affect the line emis- 
sion in the collapsing inner region). We therefore used 
the low-lying lines to constrain the other two parame- 
ters, namely Xout and M . For this we minimized the 



\ 





Fig. 2. surface as function of the water abundance in 
the outer envelope Xout and the mass accretion rate M . 
The has been obtained considering the lines with an 
upper level energy Eup lower than 142 cm^^ and for a 
central mass of 0.5 Mq . 
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obtained considering 



only water lines having a Eup lower than 142 cm~^, and 
where a is the error associated with each line flux. Fig. 
^ shows the surface as a function of Xout and M for 
Ah = 0.5 Mq . We obtained a similar plot for M* = 0.3 
Mq and the result is basically the same. As suspected, 
the chosen lines constrain relatively efficiently the two pa- 
rameters Xout and M . The minimum x^ is obtained for a 
water abundance Xout ^ 5 x 10^^ and an accretion rate 
M ~ 5 X 10^5 Mq yr-i . 

3.3.2. Central mass and water abundance in the inner 
parts of the envelope 

We then constrained the central mass M* and the abun- 
dance in the innermost parts of the envelope Xin using 
the high-lying lines. In fact to be excited, these lines re- 
quire relatively high temperatures and densities which are 
likely to be reached only in the innermost parts of the en- 
velope. Their intensities depend hence on the water abun- 
dance Xin in these parts and on the central mass M.^ . The 
X^ surface as function of these two parameters is shown 
in Fig. 1^, obtained considering the lines having an upper 
level energy Eup larger than 142 cm"-*^, and assuming M ~ 
5 X 10"^ and Xout = 5 x 10^^. In this case the x^ has a 
minimum around M* = 0.5 Mq and Xin = 5 x 10~^. 
Specifically, if we adopt a constant H2O abundance of 
5 X 10^ {Xout) across the entire envelope, the model pre- 
dicts intensities a factor between two and five lower than 
those observed (high lying lines, i.e. with Eup > 250 K). In 
other words, the observed emission can only be accounted 
for if a jump in the water abundance is introduced when 
the dust temperature exceeds 100 K, the sublimation tem- 
perature of icy mantles. This jump has to be larger than 
about a factor 10. 
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Fig. 3. surface as function of the central mass M* and 
water abundance in the innermost parts of the envelope 
Xin ■ The has been obtained considering the lines with 
an upper level energy Eup larger than 142 cm~^ and as- 
suming M = 5 X 10^^ and Xo„4 = 5 x 10^^. Note that we 
did not include the 82 and 99.5 /xm lines, which seems un- 
derestimated by our model (see text). We did not include 
the 113 /im line either, because of the blending with the 
CO J„p — 23 line, which makes the estimate of the flux 
rather uncertain. 
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Fig. 4. Ratio between the line fluxes predicted by our best 
fit model and observed ones as function of the upper level 
energy of the transition Eup . Triangles represent ortho 
H2O transitions, and squares represent para H2O transi- 
tions. Note that the model assumes an ortho to para ratio 
equal to 3. 



Parameter 


IRAS4 


IRAS16293-2422 


Mass (Mq ) 


0.5 


0.8 


Accretion rate (Mq yr^""^ ) 


5 X 10"^ 


3 X 10"^ 


Water abundance Xout 


5 X 10"^ 


3 X 10"'^ 


Water abundance Xin 


5 X lO"*^ 


3 X lO"'' 


Radius (Tdust = 30 K) (AU) 


1500 


4000 


Radius (Tdust = 100 K) (AU) 


80 


150 


Age (yr) 


1.0 X 10* 


2.7 X 10'' 



Table 2. Best fit parameters and derived quantities for 
IRAS4 and comparison with the values obtained towards 
IRAS16293-2422 bylCeccarelh et al.| (pOOOa|). 



3.3.3. The best fit model 

Assuming two identical envelopes, the best fit model is 
obtained with a central mass of 0.5 Mq , accreting at 
5 X 10"^ Mq yr-i (Table 2). 

Assuming a constant accretion rate, this gives an age 
of 10000 years, close to the dynamical age of the outflows. 
The abundance of water in the outer parts of the enve- 
lope is 5 X 10~^ and it is enhanced by a factor 10 in the 
innermost regions of the envelope, where grain mantles 
evaporate. Figure |^ shows the ratio between the observed 
and predicted line fluxes as function of the upper level en- 
ergy of the transition. The model reproduces reasonably 
well the observed water emission, with the exception of 
the lines at 99.5 /xm and 82.0 /im that seems to be un- 
derestimated (by a factor 10) by our model. Since, on the 
contrary, lines in a comparable range of Eup are well re- 
produced by our model, we think that this discrepancy is 
likely due to a rough baseline removal. Specifically, the es- 
timate of both the 99.5 and 82.0 /im line fluxes may suffer 



of an incorrect baseline removal, as the lines lie on the top 
of a strong dust feature, which covers the 80-100 /Ltm range 
(Ceccarelli, Caux, Tielens et al. in preparation). Higher 
spectral resolution observations are required to conflrm 
this explanation. Finally some unexplained discrepancy 
between the model and the observations may exist at the 
higher values of Eup . 

In the figure we also report different symbols for the or- 
tho and para water transitions respectively. In our model 
we assumed that this ratio is equal to 3. The compar- 
ison between the observations and predictions is consis- 
tent with this assumption. Plots of the predicted inten- 
sity of various lines as function of the radius are reported 
in Fig. 1^. Finally, the envelope model predicts an inten- 
sity of 1.8 X 10^^^ erg s^^ cm^^ and a linewidth of ~ 
1 Km s^^ for the ground water line at 557 GHz, equiv- 
alent to an antenna temperature of 30 mK in the SWAS 
beam (- 4'). SWAS detected a T; ~100 mK line, which 
linewidth is ~ 18 Km s~^, and self-absorbed at the rest 
velocity (Bergin et al. 2002). The observed line is undoubt- 
edly dominated by the outflow emission, and only a small 
fraction of its intensity can be attributed to the envelope 
emission, in agreement with our model predictions. 

3.4. [01], [CI I], OH and CO emission 

[01] 63 /im and [CII] 157 /tm emission is widespread, and 
probably associated with the cloud. A plausible explana- 
tion is that the two lines are emitted in the PDR resulting 
from the UV and/or X-ray illumination of this cloud from 
the several young stars that it harbors. The comparison 
of the observed fluxes with the model by Kaufman et al 
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transition, the [01] 63 /im line is relatively easily opti- 
cally thick, and an emission from an embedded source can 
be totally absorbed by the foreground material ( Poglitscli] 



et al. 



et al 



199e 



|200C) 



; paluteau et alj|1997| ; paux et al.]|l999b| ; |Vastci 



Finally, our model predicts OH and CO Jup > 14 line 
fluxes more than ten times lower than those observed. 
Note that the FIR CO hues predicted by the CHT96 
model are optically thick and not sensitive to the adopted 
abundance, and therefore increasing the CO abundance 
would not change the result. An extra heating mechanism 
is evidently responsible for the excitation of the FIR CO 
lines observed in the central position. Shocks have been 



invoked in the literature (e.g. CeccarcUi et al 



et al 
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Giannini et al. 2001), but this hypothesis has 



10 



10 



10 



10 

t3 



-12 



-13 



-14 



-15 



10 



10 



-17 




10 



100 1000 
r (AU) 



10000 



Fig. 5. Predicted intensity of various lines as function of 
the radius. The upper panel shows the water lines at 179 
/zm (solid line), 108 fim (dotted line), 75 /xm (dashed line) 
and 82 /im (dash-dotted line). The lower panel shows the 
CO Jup — 3 (solid Hue), CO J„p = 14 (dashed Hue) and 
C^^O J„p = 3 (dash-dotted line). 



( |l999| ) suggests a PDR with a density of about 10'' cm"'^ 
and a incident FUV of ~ 5 Go- This PDR would account 
for the total observed flux of the [CII] 157 /Ltm line and 
01. The parameters we derive are in agreement with those 
quoted by Molinari et al. (200C), who studied the region 
around SVS13. 

The thermal emission from the envelope predicts no 
C"*" emission, of course, as no source of ionization is con- 
sidered in the CIIT96 model. The atomic oxygen, on the 
contrary, is present all along the envelope and it is pre- 
dicted to emit 1.8 x 10^^^ erg s~^ cm^^. This is similar 
to the observed [01] 63 fim flux. The fact that we do not 
see any [01] 63 fim enhancement towards the source with 
respect to the surroundings can be explained if IRAS4 is 
well embedded in the parental cloud. Being the ground 



its own drawbacks and flaws (see Introduction). Another 
possibility is that the FIR CO lines are emitted in a su- 
perheated layer of gas at the surface of a flaring disk, as 
seen in the protostar El 29 (Ceccarelli et al. 2002), and/or 
at the inner interface of the envelope itself (Ceccarelli, 
HoUenbach, Tielens et al. in preparation). In the first case 
(disk surface) , the gas is "super-heated" because the grain 
absorptivity in t he visible exceeds the grain emissivity in 
the infrared (e.g. Chiang fc Goldreich 1997). In the latter 
case (envelope inner interface) the extra heating is pro- 
vided by the X-ray photons from the central source. A 
full discussion about the FIR CO emission origin is be- 
yond the scope of this article and we refer the interested 
reader to the above mentioned articles. 



4. Discussion 

4.1. The structure of the envelope 

The first result of our modeling is its capacity to reproduce 
the observed water emission. In the following we show that 
the derived parameter values are in agreement with pre- 
vious estimates obtained through different studies. Before 
discussing in detail this comparison we describe here the 
derived density and temperature profiles, as well as the 
gas heating and cooling mechanisms. The reconstructed 
density and temperature profiles are shown in Fig. |6[ In 
the outer region the density follows a law (static re- 
gion), while in the inner region the density is proportional 
to r^"^/^ (free-fall region), with a transition region start- 
ing at ~ 1500 AU, that covers a substantial range of radii 
and is flatter than r~^/^. The absolute density is relatively 
large, 4 x 10^ cm^^ at 1500 AU. 

The gas temperature closely follows the dust tempera- 
ture in the outer and intermediate regions of the envelope. 
At r ~ 80 AU, dust and gas decouple. This is caused by 
the evaporation of the icy mantles when the dust tem- 
perature reaches 100 K, which injects large amounts of 
water in the gas phase (about a factor ten more) increas- 
ing the cooling efficiency of the gas. The gas temperature 
drops by about 20 K and remains lower than the dust 
temperature until at r ~ 30 AU the dust FIR pumping of 
the water molecules counterbalances the water line cool- 
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Fig. 6. Density and temperature profiles of the envelope 
as computed by the best fit model. In the lower panel the 
dashed Hue refers to the dust temperature, while the solid 
line refers to the gas temperature. 



ing and couples again the gas and dust temperatures. The 
heating and cooling of the gas in the envelope follow the 
general properties discussed in CHT96 (see also Ceccarelh 
et al. 2000a). The heating is dominated by the NIR water 
pumping, gas-grains collisions and the compression of the 
gas in the innermost regions, while gas-dust collisions are 
the main heating factor in the outer parts of the envelope 
and when the dust and gas decouple at r '-^ 80 AU. The 
gas cooling is dominated by the line emission of water, 
oxygen and CO. In the outer parts of the envelope the 
cooling is dominated by the CO line emission, while in 
the intermediate region the cooling is dominated by the 
oxygen and water lines, as the CO lines become rapidly 
thick. At r < 80 AU, the cooling by the water lines takes 
over, due to the icy mantle evaporation, and dominates 
the cooling by orders of magnitude with respect to O and 
CO. Note that before the ices evaporation (i.e. at r >8G 



AU), the heating is dominated by the gas compression 
and the cooling by the II2O line emission. The increase 
of the gaseous water abundance by a factor ten causes 
an increase of about the same amount in the gas cooling 
rate, whereas the compression heating rate just increases 
by 20%. The gas heating becomes hence dominated by 
the gas-grains collisions, which tend to couple dust and 
gas. In the specific case of IRAS4, the water cooling rate 
AjjaO ^J^d gas-grains collisions heating rates r^g can be 
approximated by the following expressions, at the radius 
just before the evaporation (for the details see CHT96): 



A 



15 x(H20) /Tg 



3/2 



- 1 —3 
s cm 



(6) 



H,o = 3.4 X 10-^^ ^ ^ ; ^ erg 
5 X 10-^ V90Ky ^ 

Tdg - 3.6 X 10-16 (Td„st - Tgas) ' erg s-^cn^T^ 



The increase of a factor ten in the water cooling rate 
after the ice evaporation is only in partly counter-balanced 
by the increased heating rate due to the gas-grains colli- 
sions, and the difference between the two temperatures is 
- 20 K. 

The CO and atomic oxygen abundances are constant 
across the envelope, within the studied range, i.e. 30 to 
3000 AU. We will discuss in the next paragraph the ef- 
fect of varying the CO abundance across the envelope to 
take into account the CO depletion when the dust temper- 
ature is below the CO-rich ice evaporation temperature. 
As widely discussed previously, the water abundance un- 
dergoes a jump of about a factor ten at r < 80 AU, when 
H20-rich ices evaporate (dust temperature larger than 100 
K). 

One interesting prediction of this study is the existence 
of a hot core like region in the innermost parts of the enve- 
lope, where the dust temperature reaches the sublimation 
temperature of the grain mantles. 

4.2. CO depletion 

Here we want to address in some detail the issue of the 
CO depletion claimed in IRAS4 by comparing low-lying, 
millimeter CO transitions and dust continuum emission 
( [Blake et aLl|l995D . Blake et al. used the CO J = 3-2 
transitions of the ^^C, ^^O and ^^O isotopes and found 
that the emission is accounted for a CO abundance of 



~ 2.5 X 10~^, i.e. a factor five lower than the "canonical" 
value. Their explanation is that CO molecules freeze out 
on the grain mantles and the gas phase CO results there- 
fore depleted across the envelope. This result has been re- 
cently confirmed by JSD02, who quote CO/H2 '-^ 10^^ in 
the outer envelope of IRAS4. The question arises whether 
the CO is depleted across the entire envelope or not, as 
also remarked by JSD02. We checked with our model if 
the millimeter observations would be sensitive to higher 
abundances of CO in regions where the dust tempera- 
ture exceeds the CO-rich ice evaporation temperature, i.e. 
~ 25 - 30 K. We found that if the C^^^O J = 3-2 line is 
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optically thin, the bulk of the J — 3-2 line is emitted at 
~ 3000 AU. This is because of the excitation of the CO J 
= 3-2 line itself. In the optically thin and thermalized line 
approximation it is possible to derive the following ap- 
proximate analytical expression (strictly speaking, valid 
only for the CO millimeter lines): 



(8) 



When J = 3 the integrand has a peak around Tgas ~ 30 
K, does not vary by more than 20% with temperatures 
between 20 K and 40 K, and decreases for temperatures 
larger than ~ 40 K. As a result, the line intensity in the 
outer envelope is linearly proportional to the radius of the 
emitting gas, and therefore increases going outward, i.e. 
the emission is dominated by the outer envelope where 
CO is depleted. 

We run a case in which there is no CO depletion across 
the envelope, (i.e. the CO abundance is constant and equal 
to 10^**), a case where CO is depleted on the entire en- 
velope (i.e. the CO abundance is constant and equal to 
2.5 X 10~^), and a third case where CO is depleted only 
in the outer envelope (i.e. the CO abundance is 1 x 10^^ 
when Tdust < 30 K and 10"'' in the rest of the envelope). 
The last two cases give the same intensity (within 15%) 
on both C^^O J = 3-2 and J = 2-1 lines, while the first one 
gives an intensity larger by about a factor 3. Therefore, 
the last two cases, CO depletion across the entire envelope 
and depletion in the outer envelope only, are not distin- 
guishable by the C^^O J = 3-2 and J — 2-1 observations. 
In this sense, we think that the Jup < 3 observations can- 
not probe the innermost regions and that it is possible 
that the measured CO depletion is only relative to the 
external regions of the envelope. 

This situation would be similar to what has been 
claimed to occur in IRAS16293-2422 ( peccarelli et al.| 
2001 ), based on the indirect evidence provided by the 
D2CO emission. The D2CO molecule is considered a grain 
mantle product, as gas phase reactions seem unable to 
form an appreciable amount of this molecule. Ceccarelli 
et al. showed that in IRAS16293-2422 the D2CO emission 
originates in the region of the envelope where Tdust ^ 30 
K. The proposed interpretation is that D2CO is trapped 
in CO-rich ices that evaporate when the dust temperature 
exceeds 30 K. Hence, in IRAS16293-2422 there is an outer 
region of the envelope where CO is frozen onto the grain 
mantles (Tjust < 30 K), and a regions with T^ust > 30 K 
where CO is released into the gas phase and has the stan- 
dard ~ 10~^ abundance. A similar scenario has been also 
suggested by JSD02 for other Class sources that show 
CO depletion. 



4.3. Comparison with previous studies of IRAS4 

In this paragraph we compare our results with previous 
studies dealing in a way or in another with some of the 
issues addressed in the present study. We start with the 



recent study by JSD02, who derived the density profile of 
the IRAS4 envelopes by modeling the continuum emission 
(spectral energy distribution and 450/850 /im maps simul- 
taneously) and assuming a single power law index. They 
found a power law index equal to 1.8 and 1.3 for 4A and 
4B respectively, consistent with the Shu inside-out solu- 
tion adopted in our model. At 1000 AU they estimate a 
density equal to 6 (4A) and 2 (4B) xlO^ cm~'^, which are 
quite comparable with our estimate: 6 x 10^ cm~'^. We 
emphasize that the two methods, their and our, are to- 
tally independent, and use different data, continuum and 
line observations respectively. The fact that both predict 
approximately the same density structure in a certain way 
validate both methods, or at least increases the probabil- 
ity that both models describe reasonably well IRAS4. 

Apart from the density and temperature profiles, our 
model also constrains the water abundance profile. It is 
reasonable to ask whether our predicted water abundance 
in the innermost and outer regions of the envelope are 
realistic and if they have any support from different ob- 
servations. The situation here is somewhat complicated by 
the fact that there aren't many other independent ways 
to measure the water abundance. From a theoretical point 
of view the abundance in the outer envelope, 5 x 10^^, 
can be very well compared with chemistry model predic- 
tions (e.g. Lee et al. 1996). In this respect, the value 



that we derive is certainly not extraordinary and rather 
plausible. From an observational point of view Bergin et 
al. (2002) succeeded to detect the 557 GHz water line in 
the NGC1333 molecular cloud. They estimate the water 
abundance in the region to be ~ 10"'', with unfortunately 
a relatively large error 10) du e to the many uncertain- 
ties in the excitation of the line. Moneti et al. ( 2001 ) de- 
rived a water abundance of 3 x 10~^ in the clouds in the 
line of sight of the galactic center. These authors claim 
that this is very likely the abundance of standard molec- 
ular clouds. In summary, the water abundance that we 
find for the cold region of the IRAS4 envelope is consis- 
tent with other studies. Regarding the abundance in the 
inner region, 5 x 10~^, the value that we obtain seem to be 
lower than what expected if all the water ice is injected in 
the gas phase and a large fraction of the oxygen is locked 
in this ice. A typical water ice abundance is estimated 
around 10"'' ( fTielens et al||l99l| ). However, SWAS obser- 
vations of IRAS4 and other low mass protostars suggest 
that the water abundance in their outfiows is around 10"^ 
(Neufeld et al. 2001; Bergin et al. 2002), i.e. similar to the 
value that we find. Those estimates are very rough and 



could easily be off by a factor ten (Neufeld et al. 2000), 
as they are based on one transition only, but nevertheless 
have the advantage that the observed emission is certainly 
dominated by the outflow (the spectral resolution of these 
observations is ^ 1.2 km s~^) so there are no doubts on 
its origin. Since the water abundance in the outflow would 
be probably dominated by the grain mantles released in 
the gas phase, these observations probably measure the 
water content in the mantles, very much as our observa- 
tions measure (indirectly) the water content mantles in 
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the inner hot hke region. The two measurements seem to 
be consistent in giving a rather low value. Whether this 
validates both measures is less certain than the density 
profile case: it certainly does not discredit the two mea- 
sures. Finally, even the comparison of our estimate of the 



is smaller in IRAS4 (- 1500 AU) than in IRAS16293-2422 
(~ 4000 AU), i.e. the CO depleted part of the envelope is 
relatively larger in IRAS4 than in IRAS16293-2422. This 
may explain why the CO depletion has been observed to- 
wards IRAS 4aj^^ dvan Dishoeck 



accretion rate and central mass are in good agreement 
with the previous estimates, based on a different method 



et al 



1995; Ceccarelh et al. 2000b). 



(line profile and molecule H2CO), by Di Francesco et al 
( |2001| ). We derived M = 5 x lO^'"^ Mq yr"^ against the 
11 - 4 X 10-5 M 



i0 yr 



= 5 X 10-5 Mp 

quoted by Di Francesco et al. 



( |2001| ), and M, = 0.5 Mq against the 0.23 - 0.71 Mq . 

To conclude, these studies show that the values we de- 
rive of the four parameters of our model are plausible and 
nothing of particularly surprising, with the possible ex- 
ception of the water abundance in the innermost regions. 
In other words, if we had to choose a priori those values 
we would have chosen exactly what we found. The con- 
clusion is that it is very probable that at least most of 
the observed water emission in IRAS4 originates in the 
envelopes. If any, just a small fraction should therefore be 
associated with the outflow. Our final comment is there- 
fore that care should be taken when interpreting the ob- 
served water emission towards low mass, Class proto- 
stars as due to shocks (e.g. Ceccarelh et al. 1998, Nisini et 
al 1999, GNLOl), as we showed in two out two cases that 
the massive envelopes surrounding these sources dominate 
the water emission, just because of the large total column 
density. As a matter of fact. Class I sources, which are 
characterized by less massive envelopes, do n ot show up 
strong water emission ([Ceccarelh et alJ^OOOal) . 



4.4. Comparison with IRAS 16293-2422 

The mass and accretion rate we derived for IRAS4A and 
B are of the same order of magnitude of those found 



in IRAS16293-2422 (|CeccareUi et al.||2000aD . IRAS16293- 
2422 seems more massive (0.8 Mq ) than IRAS 4A (0.5 
Mq ), and accreting at a slightly lower accretion rate 
(3 against 5 xlQ-^ Mq yr"-'^ ). Assuming a constant 
accretion rate, those values give an age of 10000 years 
and 27000 for IRAS4 and IRAS 16293-2422 respectively. 
Hence IRAS16293-2422 seems more evolved than IRAS4. 
Moreover, IRAS4 possesses an hot core like region about 
two times smaller than IRAS16293-2422 (80 AU against 
150 AU). Ground-based H2CO and CH3OH observations 
(Blake et al. 1995; Maret et al. in preparation) confirm 
that IRAS4 is in fact colder, and therefore less bright 
in these molecular transitions than IRAS16293-2422, and 
that indeed the IRAS4 hot core like region is very small. 
This fact coupled with the larger distance of IRAS4 from 
the Sun may explain the apparent difference in the molec- 
ular emission of these two sources, which is much richer 
in IRAS16293-2422. This conclusion is also in agreement 
with the relatively higher millimeter continuum observed 
in IRAS4, which implies a larger amount of cold dust sur- 
rounding this source than IRAS16293-2422. In addition, 
the region where the dust temperature is higher than 30 K 



Finally, despite this difference in the age, the water 
abundance in the envelope is remarkably similar in the 
two sources, both in the outer part of the envelope and in 
the inner ones, where ice mantles are predicted to evapo- 
rate. This is an important piece of information, suggesting 
that the ice mantle formation in the two sources under- 
went a similar process, despite the macroscopic difference 
between the two molecular clouds which the two sources 
belong to. In the case of IRAS16293-2422, the cloud seems 
very quiescent, shielded from strong UV and/or X-ray ra- 
diation 



(.e.g. 



Castets et al. 2001) and with even evidence of 



large CO depletion ( Caux et al. 1999b). In the other case, 
IRAS4, the cloud presents caviti es excavated by th e sev- 
eral young stars of the region (e.g. Lefloch et al. 199§| ), and 
it is probably permeated by the X-rays emitted by them. 
A forthcoming study will allow to measure the II2CO and 
CH3OII abundances in the inner hot core like region of 
IRAS4 (Maret et al. in preparation) and make compar- 



isons w ith that found in IRAS16283-2422 ( [Ceccarelli et al 
^OOOb ). This study will hence help to understand in more 



detail how apparently different conditions in the parental 
clouds affect the grain mantle composition. 

5. Conclusions 

We presented a spectral survey of the protobinary system 
IRAS4 in the NGC1333 cloud, using ISO-LWS in grating 
mode. We targeted the source as well as two adjacent po- 
sitions, NE-red and SW-blue, that encompass the red and 
blue lobes of the outflow emanating from IRAS4, respec- 
tively. The three spectra are dominated by the [01] 63 /.tm 
and CII [157] /im lines, that likely originate in the PDR 
associated with the parental cloud. On the contrary, water 
emission is only detected towards the on-source position, 
whereas no significant water emission is detected towards 
the NE-red and SW-blue positions. This suggests that the 
bulk of the water emission is due to the thermal emis- 
sion of the protostellar envelopes around the two proto- 
stars. Using an accurate model of the chemistry, thermal 
balance and radiative transfer in protostellar envelopes 
(CIIT96), we modeled the water line emission due to two 
identical envelopes surrounding IRAS4A and B respec- 
tively. We found that the observations are consistent with 
the CIIT96 model, which implicitly assumes the "inside- 



out" theory (Shu 1977). The best fit of the model allows 
us to estimate the four main model parameters: the accre- 
tion rate, 5 x 10"^ Mq yr"^ , the central mass, 0.5 Mq , 
the water abundance in the outer envelope, 5 x 10^^, and 
in the inner envelope, 5 x 10"^ (this last parameter is the 
least constrained with about a factor 3 of uncertainty). 
This gives an age of 10000 years, assuming that the ac- 
cretion rate remains constant during the collapse. Based 
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on this model, we derived the density and temperature 
profiles of the gas in the envelopes. We also reviewed the 
suggestion by Blake et al. (1995) that CO is depleted by 
about a factor ten in the envelope of IRAS4. We could 
not confirm or rule out this hypothesis but caution that 
the transitions used by this study, C^^O 3-2 and 2-1, can 
hardly probe the inner regions, where the CO abundance 
may be "canonical" . 

A comparison with several previous studies of the same 
source (Blake et al (1995), Neufeld et al 2000, DiFrancesco 
et al. 2001, JSD02) shows that the derived parameters 
are reasonable and consistent with the available litera- 
ture, hence re-enforcing the thesis that the observed water 
emission is indeed due to the thermal emission from the 
envelopes. A by-product of the present study is the predic- 
tion of the existence of a hot core like region in the inner 
parts of the envelope, where grain mantles evaporate, re- 
leasing large amounts of water (about a factor ten) in the 
gas phase. Such a hot core has already been proposed to 
exist around IRAS 16293-2422, where a similar study as 
been carried out ( peccarclli et al. Z0004 |b|). Comparison 
between the two protostars, show that IRAS4 is younger 
and surrounded by a more massive envelope. This explains 
the larger continuum emission and the larger depletion 
factors observed in IRAS4. Finally, this study emphasis 
the necessity of ground based observations, where higher 
spatial and spectral resolutions are achievable. H2CO and 
CH3OH are of particular interest as they are among the 
most abundant components of grain mantles, and are 
therefore expected to evaporate in the inner parts of the 
envelope. Appropriate transitions can hence be used to 
constrain the physical and chemical condi tions in the in- 
nermo st part of protostellar envelopes (see Ceccarclli et al 
2000bD . 
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Fig. 1. ISO-LWS spectra observed towards IRAS4 on source (top line), NE-red (middle line) and SW-blue (bottom 
line). 



